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Abstract 


The phase relation between the perturbation kinetic energy (Ji ') associated with the 
tropical convection and the horizontal-mean moist available potential energy (P) associated 
with environmental conditions is investigated by an energetics analysis of a numerical 
experiment. This experiment is performed using a 2-D cloud resolving model forced by the 
TOGA-COARE derived vertical velocity. The imposed upward motion leads to a decrease 
of P directly through the associated vertical advective cooling, and to an increase of I\ f 
directly through cloud related processes, feeding the convection. The maximum K f and 
its maximum growth rate lags and leads, respectively, the maximum imposed large-scale 
upward motion by about 1-2 hours, indicating that convection is phase locked with large- 
scale forcing. The dominant life cycle of the simulated convection is about 9 hours, whereas 
the time scales of the imposed large-scale forcing are longer than the diurnal cycle. 

In the convective events, maximum growth of I\ f leads maximum decay of the per- 
turbation moist available potential energy (P f ) by about 3 hours through vertical heat 
transport by perturbation circulation, and perturbation cloud heating. Maximum decay 
of P f leads maximum decay of P by about one hour through the perturbation radiative 
processes, the horizontal-mean cloud heating, and the large-scale vertical advective cool- 
ing. Therefore, maximum gain of K f occurs about 4-5 hours before maximum decay of 
P. 
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1* Introduction 


Tropical convection occurs as a result of release of unstable energy of its embedding 
environment. The large-scale environment provides favorable thermal and moisture condi- 
tions for convection occurrence and development, on one hand, and on the other hand, it is 
adjusted by redistributing vertical thermal, moisture, and momentum structures through 
the convection. Such interaction allows us to use environmental conditions to estimate the 
properties of the convection such as the precipitation, which is a cumulus parameteriza- 
tion. Since the environmental time scales (a few days and longer) are much longer than 
the convective time scales (a few hours and shorter), the rate of production of available 
potential energy by the large-scale processes is nearly balanced by the rate of consumption 
of the available potential energy by the convection (Manabe and Strickler 1964). This 
quasi-equilibrium concept is the basic premise of the cumulus parameterization scheme 
proposed by Arakawa and Schubert (1974). The decrease of convective available potential 
energy (CAPE) that measures the thermal and moisture conditions of the environment 
often coincides with the convection development so that the CAPE and rain rate are neg- 
atively correlated (e.g., Thompson et al. 1979; Cheng and Yanai 1989; Wang and Randall 
1994; Xu and Randall 1998). The phase relation between the CAPE and rainfall may be 
related to the coupling between environmental dynamic and thermodynamic fields (Cheng 
and Yanai 1989). 

Recently, the phase difference between the CAPE and rainfall has been found as 
the evidence for relaxing the quasi-equilibrium assumption in cumulus parameterization 
(e.g., Betts and Miller, 1986; Randall and Pan 1993). The minimum CAPE occurs a few 
hours [e.g., 3 hours in GATE (Global Atmospheric Research Program Atlantic Tropical 
Experiment); 6 hours in (TOGA-COARE) Tropical Oceans Global Atmosphere Coupled 
Ocean- Atmosphere Response Experiment] after maximum rainfall. Such a phase lag is 
also demonstrated by Xu and Randall (1998) in their 2-D cloud resolving simulations. 
Xu and Randall (1998) interpreted the maximum phase lag as the adjustment time scale 
from disequilibrium to quasi-equilibrium states in the presence of time- varying large-scale 
forcing. 
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In this study, the physical processes responsible for such phase relation are examined 
through the analysis of energy conversion processes between available potential energy 
and kinetic energy in a 2-D cloud resolving simulation. We, first, establish the phase 
relation between available potential energy and kinetic energy, and use a set of energetics 
equations (section 2) to examine the essential physical processes determining the phase 
relation (section 3). The phase relation is discussed in Section 4, and the conclusion is 
given in section 5. 


2. Formulations for model, energetics, and CAPE 


2a. Model 

The cloud resolving model was originally developed by Soong and Ogura (1980), Soong 
and Tao (1980), and Tao and Simpson (1993), for studying deep convective response to 
the specified large-scale forcing. A 2-D version of the model used by Sui et ah (1998) 
and modified by Li et ah (1999) is used in this study. The governing equations with an 
anelastic approximation can be expressed as follows: 
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Here u, and w are zonal, and vertical wind components; 6 and q v are potential temperature 
and specific humidity respectively; C = (q ci q r , qg ), <? C7 9r 7 <7o and q g are the mixing 

ratios of cloud water, rain, cloud ice, snow, and graupel, respectively; p is a mean air density 
which is a function of height only; wtv is a terminal velocity which is zero for cloud water 
and ice; tt— ( p/p 0 ) K -> ft — R is the gas constant, c p is the specific heat of dry air at 
constant pressure p, and p 0 — 1000 mb ; c, e, c?, and 5 denote condensation, evaporation, 
deposition, and sublimation, respectively; Q cn — L v (c — e) + L s (d — s) + L /(/ — m) denotes 
the net latent heat release through phase changes among different cloud species, where / 
and m are fusion and melting, respectively; L 5 , and Lf are heat coefficients due to 
phase changes; Qr is the radiative heating rate due to convergence of net flux of solar and 
infrared radiative fluxes; Sc is source and sink of cloud content determined by microphyical 
processes; Ds are dissipation terms; overbar (“) denotes a zonal-mean; subscript & denotes 
an initial value, which does not vary with time; superscript 0 denotes imposed observed 
variables in the model. More discussions about the model and the model responses to 
prescribed TOGA-COARE forcing are reported in Li et al. (1999). 


2b, Energetics equations 

The perturbation kinetic energy and zonal-mean and perturbation moist available 
potential energy are respectively defined by 
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where h = c p T + L^; T is temperature; T = K0(pT) ) 1 , which is a 

parameter related to static stability; the angle bracket implies a vertical integration: 


< 0 >= r P()dz . ( 3 ) 

j ZB 

Here zq and zt are the heights of bottom and top of the model, respectively. Note that the 
horizontal-mean moist available potential energy (P) is different from the moist available 
energy (MAE) proposed by Lorenz (1978, 1979). The difference arises from the definition 
of hb. Here, hb is defined as a reference state, and is defined as the initial state (anelastic 
approximation), whereas the MAE is calculated based on the assumption that the parcel 
conserves its equivalent potential temperature in the Lagrangian frame, and equivalently, 
hb varies with time as the surface temperature and moisture vary. The vertical component 
of the moist available energy is a generalization of CAPE (Randall and Wang 1992). In 
the following derivations of energetics equations, the dissipation terms are excluded. 


An equation for the perturbation kinetic energy (K f ) can be derived by multiplying 
(2) by u f and (3) by w 1 and applying zonal-mean and the integration defined by (8) on the 
resulting equation: 
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Here, C( Jv , K f ) is the conversion between K and K f through covariance between 
perturbation zonal wind and vertical velocity under vertical shear of imposed horizontal- 
mean zonal wind, and between perturbation vertical velocities under vertical shear of 
imposed horizontal-mean vertical velocity. C{P\K ! ) is the conversion between P f and K f 
through covariance between perturbation vertical velocity and temperature. Gq v (K f ) and 
G qi (K f ) are the generation terms of K* through covariance between perturbation vertical 
velocity and specific humidity, and between perturbation vertical velocity and cloud mixing 
ratio, respectively. 


To derive the equations for the zonal-mean and perturbation moist available potential 
energy, the following equation is formed by multiplying (4) by c p n and (5) by L v and 
adding the resulting equations: 
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The equations for the zonal-mean moist available potential energy (P) and the per- 
turbation moist available potential energy (P r ) can be derived by multiplying (10) by 
c~ 1 T(h — /i&), and by c“ 1 T/i', and applying zonal-mean and the integration defined by (8) 
on the resulting equations: 

^ = C(P',P) + G r (P ) + G cn (P ) + C h (K,P) + C V (K,P), (11) 
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Here, C(P ' , P) is the conversion between P' and P through covariance between h — h b 
and convergence of vertical flux of potential temperature and moisture. Gr(P ) and G cn (P) 
are the generation terms of P through covariances between h — h b and horizontal-mean 
radiative heating, and between h — h b and horizontal- mean heating due to phase change 
of the cloud contents respectively. Ch(K,P ) and C V (K,P) are the conversion between 
K and P through covariances between h — h b and imposed horizontal temperature and 
moisture advections, and between h — h b and the horizontal-mean vertical temperature and 
moisture advections by imposed vertical velocity respectively. Gr(P') and G cn (P' ) are the 
generation terms of P' through covariances between h' and perturbation radiative heating, 
and between h' and perturbation heating due to phase changes of the cloud contents, 
respectively. G(P') is the generation term of P' . Note that C(P,P') + C(I\' ,P') + G(P') 
causes changes of P' due to the vertical advection processes. The key steps to derive (9) 
and (12) can be found in Appendix. 

2c. CAPE calculation 
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The CAPE can be calculated by 


CAPE = g [ C dpcl ^ f e ” v ^ dz. (13) 

JLFC &env{z) 

Here 6 pc i is the potential temperature of an air parcel lifted from zb to zt while not mixing 
with its environment (d env ). The air parcel is lifted dry adiabatically until it becomes 
saturated and then is lifted moist adiabatically thereafter. The level of free convection 
(LFC) is the height where 6 pc i > 6 env: z c is the level where 0 pc i = 0 env . 


The CAPE is calculated for a pseudo-adiabatic process and a reversible moist adiabatic 
process, respectively, in this study. In the pseudo-adiabatic process, an air parcel is lifted 
adiabatically while all condensed water drops out from the parcel. In the reversible moist 
adiabatic process, an air parcel is lifted adiabatically while all condensed water is kept 
in the parcel. Following Xu and Emanuel (1989), The virtual temperatures ( T vpa ) for 
the pseudo-adiabatic process and ( T vre ) for the reversible moist adiabatic process are 
respectively expressed by 
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where T p is the temperature of a pseudo-adiabatically displaced air parcel; q vs is the 
saturation specific humidity, and q w is the total water content of the air parcel. CAPE pa 
for the pseudo-adiabatic process and CAPE re for the reversible moist adiabatic process 
are calculated by using (13a) and (13b), respectively. 


3. Results 

Figure 1 shows the time evolution of vertical distribution of the large-scale vertical 
velocity and zonal wind during 19-25 December 1992 that are imposed in the model. 
Strong upward motions with maxima of 15 — 25 mb hour occur on late 20 December, 
and during the early mornings of 23 and 25 December between 400 and 500 mb. The 
latter two maxima are quasi-two-day oscillations (Takayabu et al. 1996) in the convective 
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phase of an intraseasonal oscillation during COARE. Two less intense upward motion 
centers appear during the nights of 19 and 21 December. The occurrence of maximum 
upward motion each night is consistent with the diurnal signals observed by Sui et al. 
(1997). The large scale zonal wind in the lower troposphere (below 700 mb) are westerly 
that strengthens to 10 ms” 1 around 23 December. The mid-troposphere has an easterly- 
westerly wind oscillation with maximum easterly wind of -10 ms~ l at 500 mb on 20 
December. The upper troposphere (above 250 mb) is dominated by easterly winds. As 
mentioned previously, the model is also forced by the observed horizontal temperature and 
moisture advections (not shown), which have much smaller amplitudes than the vertical 
advections respectively. 

Figure 2a shows lag correlation coefficients between zonal-mean CAPE and rain rate. 
Positive lag hour denotes that CAPE leads rain rate. The maximum lag correlation coeffi- 
cients between zonal-mean CAPE and rain rate indicate that the CAPE reaches maximum 
about 3-4 hours before the maximum rain rate. The minimum lag correlation coefficients 
indicate that the CAPE reaches minimum about 2 hour after the maximum rainfall. Both 
maximum and minimum are above 99% confidence level. The phase difference between 
maximum and minimum correlation coefficients is about 5 hours. Since a significant spec- 
tral peak appears at 9 hours by the power spectrum analysis of the hourly rain rate (not 
shown), the phase difference is about the half of the time scale of the simulated convec- 
tion. Figure 2a also show that the lag correlation coefficients for C APE re and CAPE pa 
are similar. 

Since the model is forced by imposed vertical velocity, the relationship between energy 
and imposed vertical velocity is first analyzed. Figure 2b shows lag correlation coefficients 
between P and < w° > (solid line), and between K f and < w° > (dashed line). Statisti- 
cally significant lag correlation coefficients display that maximum K f lags imposed upward 
motion (positive < w° >) by 1-2 hours whereas minimum P lags upward motion by about 
6 hours. This suggests that the I\ r leads P by about 4-5 hours, which is about the half 
of the time scale of the simulated convection. The statistically significant relationship can 
be also shown by lag correlation coefficients in Figure 3a that the imposed upward motion 
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leads minimum (maximum decrease of horizontal-mean moist available potential en- 
ergy) by 3 hours, whereas it lags maximum (maximum increase of perturbation kinetic 
energy) by 1-2 hours. Thus, minimum lags maximum by 4-5 hours. The imposed 
large-scale upward motion with maximum A^- and K' indicates that convection is phase 
locked with the imposed large-scale upward motion. The negative lag correlation coeffi- 
cient between and < w° > in Figure 3 means that the imposed large-scale downward 
motion leads maximum by about 3 hours. Thus, the imposed downward motion results 
in a build up of P, and provides a favorable environmental conditions for occurrence of 
convection. 

The phase relation between A^- and A^- is also linked by local change of perturbation 
moist available potential energy Ajp. Minimum A£ lags minimum Ajj- by about 1 hour, 
and minimum Aj^- lags maximum by about 3 hours (Fig. 3b), so that minimum Ap- 
ings maximum A^p- by about 4 hours. This is a statistically significant phase relation 
consistent to that shown in Figure 2b, although there are other two lag correlation coef- 
ficients between A^- and A^- that are above the 99 % confidence level. Two minimum 
lag correlation coefficients are 9-hour apart, indicative of the dominant life cycle of the 
model convective events. The maximum lag correlation coefficient appears between the 
two minimum lag correlation coefficients, indicating that A^p- reaches maximum about one 

nr'/ 

hour before maximum g jrf- 

To further examine the dominant physical processes determining the phase relations, 
the lag correlation between each term of A£ [eq. (11)] and A^p- , and between A^ and each 
term of A^- [eq. (12)], and the lag correlation between each term of A^- [eq. (12)] and 
A^A-, and between A^- and each term of A^- [eq. (9)] are respectively plotted in Figures 
4 and 5. Figure 4a shows that only the zero-hour lag correlation coefficient between 
C V (K,P) and A^- is marginally around the 99 % confidence level, and above the 95 % 
confidence level. The term C V (K,P) of Aj^ is a major component that contributes to the 
maximum positive zero-hour lag correlation between A£ and A^-. The term C V (K,P ) is 
contributed by the vertical temperature and moisture advections respectively [see (lie)]. 
Further analysis shows that the lag correlation coefficient between vertical temperature 
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advection and imposed vertical velocity has the same sign with those between sum of 
vertical temperature and moisture advections and imposed vertical velocity, whereas the 
lag correlation coefficient between vertical moisture advection and imposed vertical velocity 
has the opposite sign (not shown). This indicates that vertical temperature advection 
determines the conversion term C V (K , P ). The imposed upward (downward) motion causes 
the vertically advective cooling (warming), which results in the loss (gain) of horizontal- 
mean moist available potential energy through the conversion term C V (K,P). Figure 
4a. also shows that maximum zero-hour lag correlation coefficient between G cn (P ) and 

is slightly less than that between C V (K,P) and The conversion G cn (P ) carries 
convective signals as shown in Figure 4a that the maximum lag correlation coefficient is 
about 9-hour apart.. 

Figure 4b shows that the lag correlation coefficients between C(P,P') and and 
between Gr(P') and and between G(P') and are above 99 % confidence level. The 
terms C(P, P ') and G(P') have the same orders of magnitude (not shown), but they have 
the opposite signs (Fig. 4b) so that they cancel each other in large part. In addition, the 
amplitude of the term C(K',P') is smaller than those of the terms C(P,P') and G(P'). 
As a result, the lag correlation coefficient between C(P,P') + G{P ' ) + C(K',P') and 
becomes statistically insignificant. This suggests that the vertical perturbation advection 
processes do not play important roles in determining phase of Therefore, the term 

Gr(P') makes that leads by about 1 hour. The radiative cooling with positive 
h 1 and radiative warming with negative h' cause the decrease of perturbation available 
potential energy through the conversion term Gr(P'). 

Figure 5a shows that three maximum negative lag correlation coefficients and one 
maximum positive lag correlation coefficient are above 99 % confidence level. Again, small 
amplitude of C(K',P r ) in contribution to and cancelation between C(P,P') and 

G(P') make the vertical perturbation advection processes less important in determining 
the phase of The term G cn (P') plays a crucial role in controlling the phase of as 
shown in Fig. 5a where the lag correlation coefficient between G cn {P ' ) and has the 
maximum negative value at -2 hour, and is statistically significant. This suggests that the 
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maximum loss of perturbation moist available potential energy and the maximum gain of 
perturbation kinetic energy are linked by the term G C n(P r )- Since L s > the heating 
released by deposition with positive h! causes the loss of perturbation moist available 
potential energy. Maximum G cn (P r ) occurs about 3 hours after the maximum It is 
important to notice that the maximum K f also occurs about 3 hours after the maximum 
(Figs. 2b and 3a), indicating the minimum G C n{P r ) coincides with the strongest 
convection. Thus, the 3-hour of phase difference between and is the time length 
for convection developing to the strongest strength. 

Figure 5b shows that the lag correlation correlation coefficients between C(P\K f ) 
and and between G qv {K f ) and are similar, and their maximum negative values at 
-3 hours are barely above 99 % confidence level, which contribute to maximum negative 
lag correlation coefficient between and at -3 hours (Fig. 3b). The terms G q i(K f ) 
and C(J\ , K') have the opposite signs with the terms C(P\K f ) and G qv (K f ) (Fig. 5b), 
the magnitude of C(P f y K') is larger than the other three terms (not shown) so that 

r\ t r/ 

the term C(P\K r ) determines Covariance between the temperature and vertical 

velocity perturbations determines the local change of the perturbation kinetic energy. The 
thermally direct circulation of the upward motion with the higher temperature and the 
downward motion with the lower temperature converts the perturbation moist available 
potential energy to the perturbation kinetic energy, feeding the convection. 


4. Discussion 

Figure 6 summaries the phase relations between the convection and its embedding 
environment. The imposed large-scale downward motion yields a growth of P by the asso- 
ciated vertical advective warming [ C V (I\ J P ) > 0], building the favorable environment for 
occurrence of convection. The near simultaneous occurrence of maximum K , and 

large-scale upward motion, implying that convection is phase locked with the large-scale 
forcing. The life cycle of the simulated convective events (about 9 hours) is much shorter 
than the time scales of imposed large-scale forcing (longer than the diurnal cycle). In the 
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convective events, maximum leads maximum by about 3 hours through perturba- 
tion cloud heating [G cn (P')] and the vertical heat transport by perturbation circulations 
[C(P\K ! )\. Maximum also leads maximum K f by about 3 hours, indicating that the 
3-hour is the time required by convection to reach the maximum strength. Minimum 
leads minimum by about one hour through perturbation radiative processes [Gr(P')] 
and the horizontal-mean cloud heating [G cn (P)], and the large-scale vertical advective 
cooling. Consequently, maximum leads minimum ^ by 4-5 hours (about the half of 
the convective life scale). 

The phase difference between perturbation kinetic energy associated with convection 
and its embedding environment associated with horizontal-mean moist available poten- 
tial energy indicates that the generation of environmental unstable energy by large-scale 
processes is not simultaneously balanced by its destruction by convection. The minimum 
horizontal-mean moist available potential energy occurs 4-5 hours after maximum pertur- 
bation kinetic energy, and the phase lag is about the half of the convective time scale. This 
provides a concrete evidence for the adjustment from disequilibrium to quasi-equilibrium 
states proposed by Xu and Randall (1998). The results shows the convective time scale is 
related to the cloud microphysical processes, convective-radiative interactions, dynamic- 
thermodynamic coupling inside the convective system, suggesting that convective time 
scale (as well as the phase lag) may depend on characteristics of convection (cloud type). 
When the phase lag is included in cumulus parameterization, the convective time scale 
is introduced in the large-scale environment in the general circulation simulations. The 
scale interaction may have accumulated effects on the large-scale variabilities ranging from 
diurnal, intraseasonal, interannual time scales. This needs careful consideration in cumu- 
lus parameterization scheme, and more sensitivity experiments by the general circulation 
model. 

5. Conclusion 

Energetics analysis is carried out with a 2-D cloud resolving simulation to deter- 
mine the physical processes responsible for the phase difference between convection and 
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its embedding environment. The cloud resolving model is forced by imposed time- varying 
horizontal-mean vertical velocity as well as zonal wind derived from the TOGA-COARE 
data set for a selective six day period. The imposed vertical velocity serves as external 
forcing in this particular model setup. Lag correlation analysis shows that the maximum 
perturbation kinetic energy associated with the simulated convective events and its max- 
imum growth rate lags and leads the maximum imposed large-scale upward motion by 
about 1-2 hours respectively, indicating that the convection is phase locked with the im- 
posed large-scale forcing. The imposed large-scale vertical velocity has the time scales 
longer than the diurnal cycle, whereas the simulated convective events have the dominant 
time scale of about 9 hours. The imposed large-scale upward motion leads to decrease 
the horizontal-mean moist available potential energy by the associated vertical advective 
cooling, providing the favorable environment for convection development. 

With the maximum heating released by phase changes of cloud contents associated 
with the strongest convective events and the maximum vertical heat transport by pertur- 
bation circulations, maximum growth of perturbation kinetic energy leads maximum loss 
of perturbation available potential energy by about 3 hours. With the maximum vertical 
advective cooling associated with upw r ard motion, the maximum horizontal-mean cloud 
related heating, and maximum perturbation radiative processes, maximum loss of pertur- 
bation moist available potential energy leads maximum loss of the horizontal-mean moist 
available potential energy by about one hour. Consequently, the maximum gain of pertur- 
bation kinetic energy leads the maximum loss of horizontal-mean moist available potential 
energy by about 4-5 hours (about the half of the time scale of life cycle of the simulated 
convection). 
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APPENDIX 


Relations used in derivation of energetics equations 
The following relations are derived to obtain (9): 
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where mass continuity and zero vertical velocity at top and bottom of the model atmosphere 
are applied. The following relations are derived to obtain (12): 
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Figure Captions 

Fig. 1 Time evolution of (a) vertical velocity (mb hour ~ l ), and (b) zonal wind (m 3 " 1 ) taken 
from the TOGA COARE (Sui et al. 1997) for a six-day period. Downward motion in (a) 
and westerly wind in (b) are shaded. 

Fig. 2 Lag correlation coefficients between horizontal-mean CAPE and rain rate in (a), 
and between P and < w° > (solid) and between I\ f and < w° > (dashed) in (b). Solid 
and dashed lines in (a) denote cases for a pseudo-adiabatic process and a reversible moist 
adiabatic process respectively. Lag correlation coefficient curves above upper light dot line 
or below lower light dot line exceed 99 % confidence level. 

Fig. 3 Lag correlation coefficients between and < w° > (solid) and between and 
< w° > (dashed) in (a), and between 4™ and (solid) and between 4^- and 
(dashed) in (b). Lag correlation coefficient curves above upper light dot line or below 
lower light dot line exceed 99 % confidence level. 

Fig. 4 Lag correlation coefficients between each term of ^ [Eq. (11)] and in (a), and 
between ^ and each term of [Eq. (12)] in (b). Lag correlation coefficient curves 
above upper light dot line or below lower light dot line exceed 99 % confidence level. 

Fig. 5 Lag correlation coefficients between each term of [Eq. (12)] and in (a), and 
between and each term of [Eq. (9)] in (b). Lag correlation coefficient curves above 
upper light dot line or below lower light dot line exceed 99 % confidence level. 

Fig. 6 Schematic diagram for the summary of phase relationship between the convection and 

r\ n r\ ry/ rj T .-'f 

its embedding environment in term of and 
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